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Abstract. - We develop a theory of effects of electron-electron collisions on the Dyakonov-Perel' 
spin relaxation in multi- valley quantum wells. It is shown that the electron-electron scattering rate 
which governs the spin relaxation is different from that in a single-valley system. The theory is 
applied to Si/SiGe (OOl)-grown quantum wells where two valleys are simultaneously populated by 
free carriers. The dependences of the spin relaxation rate on temperature, electron concentration 
and valley-orbit splitting are calculated and discussed. We demonstrate that in a wide range of 
temperatures the electron-electron collisions can govern spin relaxation in high-quality Si/SiGe 
quantum wells. 
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Introduction. — Electron spin dynamics is among 
the most rapidly developing branches of the modern solid 
state physics due to the rise of spintronics [1,2]. The 
prospects of spintronics which aims at the utilization of 
electron spin on equal grounds with its charge in novel 
Semiconductor devices are related with the possibilities to 
create, control and manipulate the electron spins. The 
Understanding of microscopic mechanisms of electron spin 
decoherence and relaxation is, hence, of high importance. 
' The main mechanism of electron spin relaxation in bulk 
semiconductors and semiconductor quantum wells (QWs) 
is Dyakonov-Perel' (or precession) mechanism [3,4]. It 
is connected with the spin-orbit splitting of the conduc- 
tion band states which acts as a wavevector (fe) depen- 
dent effective magnetic field with the Larmor precession 
frequency fJfc. Such an effective field arises only in non- 
centrosymmetric systems, the most widespread examples 
of them being bulk III-V semiconductors and QWs on 
their base. Although bulk Si and Ge crystals possess 
an inversion center, it has been demonstrated experimen- 
tally [5, 6] that the one-side modulation-doped Si/SiGe 
QW structures exhibit the Rashba effect and, in these 
structures, the electron spin relaxation is governed by pre- 
cession mechanism as well. Recently, a theoretical estima- 
tion for the electron spin-orbit splitting in Si/SiGe het- 
erostructures have been obtained by using the empirical 
tight-binding model computation [7,8]. 

The electron spin precession in the effective magnetic 



field is interrupted by the scattering events which change 
randomly the electron wavevector and, hence, the direc- 
tion of the spin precession axis. Thus, the spin relax- 
ation rate can be estimated as (f^fcT) where angular 
brackets denote the averaging over the electron ensem- 
ble and r is the microscopic scattering time. Hence, the 
spin relaxation is slowed down by the scattering. It is 
evident that any momentum scattering process such as 
interaction of an electron with static impurities, interface 
imperfections or phonons stabilizes the spin. It is much 
less obvious that the electron-electron scattering can also 
suppress the Dyakonov-Perel' spin relaxation contributing 
additively to [9-13] and making the time t different 
from the momentum relaxation time. Indeed, it does not 
matter whether the electron wavevector is changed in the 
process of momentum scattering, due to the cyclotron mo- 
tion or as a result of collision with other electrons [9] . It is 
established that nothing but an inclusion of the electron- 
electron scattering allows one to describe the temperature 
dependence of spin relaxation rates in high-quality GaAs 
QWs [12]. 

Here we address the electron-electron scattering effects 
on spin relaxation in Si/SiGe quantum wells. Their spe- 
cific feature is the presence of several valleys [two in case of 
(OOl)-grown QWs] populated by electrons. The Coulomb 
scattering cannot transfer an electron from one valley into 
another although electrons from different valleys can in- 
teract with each other. We show here that the micro- 
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scopic scattering time t determined by electron-electron 
collisions in the multi-valley band system is different as 
compared with the single-valley case studied previously. 
The difference is related not only to the non-equal Fermi 
energies in the single- valley and multi-valley systems with 
equal electron densities but also to the different screening 
of Coulomb interaction in single- and multi-valley bands. 

Model. — To be specific we consider Si/SiGc QWs 
grown along the axis z || [001]. The conduction band 
states arc formed from electron states in two A val- 
leys with the extrcma iiiTo = (0, 0, ±i^o)i where Kq w 
0.8 X 2'K/aQ and oq is the lattice constant. The electron 
reflection from the QW interfaces is accompanied by the 
intervalley transfers —K^ Kq and vice versa which re- 
sults in the valley-orbit splitting and formation of two sub- 
bands j = ±, the lower subband j = — and the higher one 
j ~ +■ The valley-orbit splitting Avo, depends on the 
QW width and interface properties. It may reach several 
mcV in relatively thin quantum wells [7, 14]. The electron 
eigenstates |fc, j) arc superpositions of single-valley states 
and, in the envelope-function approach, can be written as 

vl/^.(r) = e^^'^-^+'^yy^CM^)[4ly^Ko + c^J^k.] . (1) 

Here 4'±Ko the scalar bulk Bloch functions at the 
two extremum points ±Kq, kx,ky are components of the 
two-dimensional wave vector fc _L z, (f{z) is the single- 
valley envelope function calculated neglecting the interval- 
ley mixing and the spin-orbit interaction, cj^jf^ are coordi- 
nate independent scalar coefficients, |c^|, |^ + [c^^fj, |^ = 1, 
and Cs is a constant spinor describing the electron spin 
state. In QWs with asymmetric heteropotential (or with 
odd number of Si monoatomic planes) each of the sub- 
bands is split with respect to electron spin. The typi- 
cal values of the spin-splitting have /ieV range, i.e., they 
arc much smaller than the valley-orbit splitting. Conse- 
quently, the electron Hamiltonian is decomposed into two 
partial spin-dependent Hamiltonians 



2m* 



± 



2 



(2) 



describing electrons in each of the valley-orbit-split sub- 
bands. Here cr is the vector composed of Pauli matrices 
and Jlfc is the angular frequency describing the spin split- 
ting. The comparison of theoretical estimations and ex- 
perimental data [8,15,16] shows that in the state-of-the-art 
samples the spin splitting is isotropic in the QW plane and 
has a symmetry of the Rashba type, il^''-' = Pjiky, —kx, 0) 
and O^"''' = |r2^"'''| = \l3j\k. The arrangement of electron 
states is schematically shown in Fig. 1. 

The kinetic theory of spin relaxation in Si/SiCe QWs 
is developed within the density matrix method. It is as- 
sumed that the valley-orbit splitting Avo can be compa- 
rable with characteristic energy of electrons and exceeds 
by far the inverse scattering time. In this case elements 
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Fig. 1: Schematic subband structure in an n-doped Si/SiGe 
QW. The valley-orbit splitting, Avo, and spin splitting, l3+k 
and P-k, are shown not to scale. Inset illustrates population 
of the subbands j = ± by electrons, fj,+ and fx- are the chemical 
potentials referred to the subband bottoms. 



of the density matrix nondiagonal in the subband indices 
j 7^ j' can be disregarded whereas no restrictions are im- 
posed on the density matrix in the spin subspace. Within 
each subband the spin-density matrix can be recast as 



Pk 



■Ik 



U = ±) 



(3) 



where /^"'•' is the average occupation of the k state in the 
subband j , s^-'^ is the average spin in this state, the symbol 
of the unity 2x2 matrix is omitted. 

The kinetic equation for the spin density matrix can be 
represented as a set of equations for the scalar fj^^ and 
pseudovector 

4'^ as follows 



dt 



ds 



dt 



(4) 
(5) 
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Here rjp'j^ is the effective field arising from the Hartrec- 
Fock interaction in the spin-polarized electron gas [10, 17]. 
The scalar and vector electron-electron collision inte- 
grals, intra-valley (Qjf^l/, s}, Q^''''{s, /}) and inter-vallcy 

(Q'-'^I/, s}, Qfc '{s, /}), are described in the next section. 

Intra- and inter-valley interaction. The colli- 
sion integrals in Eqs. (4) and (5) describe the electron- 
electron scattering processes 



(6) 



where Si^s'^ etc. are the electron spin components ±1/2. 
Because of a long-range character of the Coulomb in- 
teraction Vc, the intervalley scattering accompanied by 
transfer of the wavevector ^ 2Ko is strongly suppressed, 
and one can exclude from consideration any contributions 



due to the matrix elements (fc^, k'y, - 
{k'^,k'y,Ko\Vc\k^,ky,-K„). 
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Fig. 2: Illustration of the direct (a) and exchange (b) Coulomb 
scattering between particle 1 with the spin s in the subband j 
and particle 2 with the spin s' and in the subband j' . 



As a result, the effective matrix element describing the 
process (6) with allowance for the indistinguishability of 
the carriers reads (c.f. [10,18]): 



M{j2ps2; j2P' S2\jiksi; j[k' s[) = 6k+k' ,p+p' 



(7) 
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where Vk-p is the Fourier-transform component of the 
quasi- two-dimensional Coulomb potential. Figures 2(a) 
and 2(b) illustrate two contributions to the scattering pro- 
cess with a pair of electrons in the final state with the wave 
vectors p and p'. The first term in Eq. (7) is a Coulomb 
interaction where the first electron changes its wave vector 
from k to p while the second electron exhibits the scatter- 
ing k' p' . The second term results from the scattering 
k ^ p' and k' p, it enters Eq. (7) with the opposite 
sign. In the classical physics, the total effective cross- 
section is proportional to the sum |Vfe_pp + |Vfe_p'p [19]. 
In quantum mechanics, for two electrons which have the 
same spin, Si = S2, and occupy the same subband, ji = j2, 
the cross-section has an additional interference term pro- 
portional to Vk-pVk-p' [20]. Note, that a simple form of 
the above equation stems from neglecting the spin-orbit 
interaction in the processes of scattering [21]. 

The collision integrals in the kinetic equations are de- 
rived by using the standard diagram technique [10] and 
Eq. (7). Here we consider the experimentally typical sit- 
uation of weak spin polarization, <ti fk^ (although 
in GaAs the realization of a remarkable optical orienta- 
tion of electron spins is also possible [17]). In this case 
the Hartree-Fock terms ^IPk kinetic equations (5) 
are unimportant and can be neglected. Let us present the 
collision integrals Qk\f,s} and Qk''{f,s} in Eq. (4) in 
a convenient form 



(3) 
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respectively. Here Ej: is the spin-independent part of the 



electron energy equal to h^k^ /2m* ± Avo/2. The above- 
defined scalar functions take the form 



Plf,I'pp' = i2V^%-Vk-p>Vk-p) 



(8) 



for the intra-subband scattering, and 



P, 



kk'pp' 



(9) 



for the subband-subband scattering, similarly to the case 
of electron-hole scattering and electron-electron scattering 
in a quantum well with several occupied size-quantized 
subbands [23] . It is worth mentioning that for the scatter- 
ing between different particles (e.g. electrons and ions in 
plasma) the scattering rates are by the factor of 2 smaller 
as compared with those given by Eq. (9) because of the 
absence of the contribution given by Fig. 2(b). 

For the pseudovector collision integrals Q'--'''{s,/} and 
Qfc^ {*'/}' similarly introduce the pseudovectors 
Pkk'vv' '^i^d Pkk'r,r,> which are given, respectively, by 



c(3) 



(2^fc% - Vk-pVk-p')[s2^F,{k';p,p') - s'jl>F,ip';k,k' 



-14-pVfc_p.[s^?F,(fe;p,p') - s'^>F,{p';k,k')] 



(3)- 



2V^-P s^^F,{k'-p,p') - 4^)F,(p'; fc, k') 



(10) 



where F.^k.-MM) = /^^(l - - /^^) + /^f /^f and 
F,{k,, k2,k,) = fi-^\l - /(7') + 7(7^1 Simi- 

lar collision integrals for subband-subband scattering were 
derived in Ref. [24] for the electron-electron collisions in 
GaAs quantum well with F and L occupied valleys. 

Before turning to the spin relaxation times we discuss 
the screening of Coulomb potential in a multivalley sys- 
tem. Assuming that the QW width is small enough to per- 
mit the electrons to be treated as strictly two-dimensional, 
the Fourier transform of Coulomb potential may be writ- 
ten approximately as, e.g., Refs. [23,25], 



Va 



27re2 



S ie{q + Qs) 



(11) 



where e is the elementary charge, S is the normalization 
area, a: is the static dielectric constant, and Qs is the in- 
verse screening length given by 



2m*e2 



3 



1 



(12) 



Here the summation is carried out over occupied sub- 
bands, fcs is Boltzmann's constant, T is the absolute tem- 
perature, is the chemical potential of electrons referred 
to the bottom of the j-th subband, see inset in Fig. 1. In 
the limit of non-degenerate electrons, exp (— /ij/ZcsT) 3> 
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1, and the screening is negUgible. If electrons are strongly 
degenerate, exp (— /ij/fc^T) <^ 1, each occupied subband 
yields the same contribution 2m* / (ech^) and the total 
inverse screening length increases proportionally to the 
number of occupied subbands. 

Spin relaxation times. — Kinetic equations (4), (5) 
are solved following the standard procedure [10]. We con- 
sider the equilibrium electron distribution with /j^"'' = 
{exp [(E'^-'^ — jj,j)/kBT] + 1}~^ and seek the spin distri- 



bution function s^^^ in the form 



6s 



U) 



(13) 



Here is a quasi-equilibrium axially-symmetric spin dis- 
tribution function related to the initially created total elec- 
tron spin in the j-th subband by S^-'' = J2k ^^"^ '^'^fc ' 
is a non-equilibrium correction resulting from the electron 
spin precession around the vector fl^j^"^ . Below we assume 

ri^-'V <C 1 (the collision dominated regime) where r is 
the typical scattering time. This condition is surely satis- 
fied in Si/SiGe QWs [5,6,15]. Since the collision integrals 



Qk^ {Ss, f} and Q)^' {5s, f} conserve the angular depen- 
dence of Ss^^^ one can present this correction as follows 



,0) 



6s 



where F^"'-' is a function of fc = ]fc]. It can be found from 
the solution of linearized Eq. (5). For the Rashba-like spin 
splitting we eventually arrive at 



1 



' s,zz 



k 



n'^' F, 



(14) 



and Ts^xx — 

ation time in the j-th subband for the spin oriented along 
the a axis. 

In the limits of degenerate and non-degenerate statistics 
it is instructive to introduce an effective scattering time 
T* in the jth subband defined by 

^ - ^ (15) 



's,yy 



'iTs Jz, where Ts%a is the spin relax- 



' s.zz 



3 3 



where the characteristic spin precession frequency Qj = 

Pjk'^P for a degenerate electron gas and fij = PjkT for a 

non-degenerate gas, kp is the Fermi wavevector at zero 
temperature in a given subband, and fc^ is the thermal 
wavevector \f2/rrtFkBT /fi. In fact, the time t* is a mi- 
croscopic electron-electron scattering time governing the 
Dyakonov-Perel' spin relaxation in each subband. Com- 
paring Eqs. (14) and (15) we obtain 



.0) 



2^k 



(3) 



(degenerate electrons) , (16) 



1.2 k 

k 



(3) 



(non-dcgencrate electrons) . (17) 



Results and discussion. — Below we present ana- 
lytical and numerical results for the microscopic scatter- 
ing times r* which govern Dyakonov-Pcrcl' spin relaxation 
in multivalley QWs. In order to emphasize the role of 
electron-electron interaction the effects of single-particle 
momentum scattering are ignored, they can be taken into 
account by inclusion into the right-hand side of kinetic 
equation (5) the collision term —6sk/Tp, where Tp is the 
momentum scattering time. 

For the non- degenerate electron gas, one can neglect 
the screening of the electron-electron interaction and the 
collision integrals describing intra-subband and subband- 
subband Qfe{s, /} and Qk{s, /} differ, apart from the 
terms Vfe_pVfc_p', by a common factor, fki^/ffJ''\ re- 
sulting from different populations of the valley-orbit split 
subbands. The inverse microscopic scattering time r*^^ 
has two additive contributions caused by the collision of 
electrons within the same subband and electrons in dif- 
ferent subbands each of those being proportional to the 
number of electrons in a given subband. Neglecting the 
terms Vk~pVk-p' and making use of the results for a single 
valley [9, 10] we have 



* _ * _ (B) 



(18) 



where Tee is the electron-electron scattering time which 
governs spin relaxation in the single valley structure occu- 
pied by electrons with total concentration N = + N- , 



_ hx^ksT 



(19) 



and / is a numerical factor which, for strictly two- 
dimensional electrons, equals to « 0.027 [9, 10]. The 
scattering times in the subbands are, therefore, the same, 
since it does not matter whether an electron scatters by 
an electron in the same or in the other subband. The spin 
relaxation times in the valley-orbit split subbands are dif- 
ferent only due to the difference of the spin splittings in 
the subbands. Note, that the allowance for the interfer- 
ence contributions Vk-pVk-p' results in a shght (~ 4%) 
increase of the constant / for the intrasubband interac- 
tion [9, 10], hence, these terms can be safely neglected. 

Now we turn to low temperatures where the electrons 
are degenerate. Figure 3 depicts the dependence of the 
scattering times r|. on Avo related to the Fermi energy 
Ep of electrons of the same concentration populating a 
single valley. In this case the electron-electron collisions 
are suppressed due to the Pauli principle and. moreover, 
the screening parameter qs is not negligible. This gives 
rise to two additional competing factors which have effect 
on the difference between the scattering rates in singlc- 
and two-valley systems. First, the electrons are redis- 
tributed between valleys which results in a decrease of 
electron concentration in each valley and, consequently, in 
an enhancement of the scattering rate due to reducing the 
Pauli blocking. Second, the screening efficiency increases 
and, therefore, the scattering rates are decreased. Due to 
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Fig. 3; Electron-electron scattering times t± in a two-valley 
QW as a function of the valley-orbit splitting. The times are 
presented in units of the similar scattering time in a single val- 
ley with the same carrier density and temperature. Black solid 
line corresponds to the lower, more populated, valley while red 
dashed line describes the upper, less populated, valley. The 
calculation is carried out for the degenerate electron gas, tem- 
perature T = 8.2 K, iV = 2 X 10^^ cm"^, Ef is the Fermi energy 
in the single valley with the same concentration A'^. Other pa- 
rameters used in the calculation correspond to Si/SiGe QWs: 
ae — 12, and m* = 0.191mo, where mo is the free electron 
mass. 



the competition between these two factors the electron- 
electron scattering time can be both longer and shorter in 
a two- valley system as compared with a single valley. 

A simple analytical expression for the electron-electron 
scattering time can be derived in the absence of intervalley 
mixing in which case iV+ = N_ = N/2, fx^ = n_ = Ep/2 
and = = T*. Let us introduce the electron-electron 
scattering rate governing the spin relaxation in a single- 
valley system with the degenerate electrons of the total 
density JV [10] 



1 



' ee ^ 

In a two- valley system at, fc^T <C /x, one has 



(20) 



(21) 



The factor 1/4 in Eq. (21) results from allowance for the 
valley-valley interaction and takes also into account that 
the Fermi energy in each valley is twice smaller as com- 
pared with the single-valley system with the same to- 
tal density. The factor J describes the modification due 
to allowance for the screening. In the limit of kp = 

\J m*Epl'fi' <?C (7s, i.e., where the screening is so strong 
that the electron-electron interaction is effectively short- 
range, J = 4 since the inverse screening length is twice 
smaller as compared with single valley system and, hence, 
the scattering probability decreases by a factor of 4. In 
real QWs kp and qs can be comparable [26] and J ranges 
from 1 to 4 depending on the electron concentration. For 



Fig. 4: Electron-electron scattering times as a function of tem- 
perature calculated from Eq. (16) for single valley (curve 1) and 
two-valley (curve 2) quantum wells. The valley-orbit splitting 
is set to zero. The electron concentration N = 2 x 10^^ cm^. 
Other parameters are the same as in caption to Fig. 3. 



the parameters used in the calculation of Fig. 3 the factor 

/ p'. 

J « 2.2 and, at Avo <C Ep, the ratio t* free is close to 
0.55. In QW structures with a large number of unmixed 
valleys, n„ > 1, like (lll)-grown Si MOSFET structures 
and degenerate electrons, the scattering time t* increases 
oc riy due to the competing effects of enhancing screen- 
ing and decreasing Pauli blocking. With the increasing 
valley-orbit splitting, the electron-electron scattering time 
in the lower valley, t_ , becomes longer and the scattering 
time r_|_ shortens. This is a result of electron redistribu- 
tion downward to the lower subband and an enhancement 
of Pauli blocking there. In the upper subband the electron 
density decreases and the Pauli blocking becomes weaker. 
If all the electrons fill the lower subband the scattering 
time, T_, rapidly drops because the screening parameter 
Qs reduces by a factor of 2 and approaches the single- valley 
value. Fig. 3. One can also see from this figure that the 
electron-electron scattering time can be both shorter 
and longer than that for the single-valley system. 

For non-zero valley-orbit splitting, the spin relaxation 
times of the electrons in the two subbands j = ± can be 
different due to the following reasons: (i) difference of the 
electron-electron scattering times 7^ t1, (ii) difference 

of the Fermi wavevectors /c^-* and (iii) difference of the 
spin-splitting constants /?+ f3-. Weak intervalley scat- 
tering characterized by the time t„ ^ , Tp may lead to 
the efficient intermixing of spins in different valleys. The 
observed spin relaxation time for the spin along one of the 
main axes a is, hence. 



T — 



2t(+) r^-) 

' s,aa ' s.aa 
I s.aa \ ' s.aa 



.(±) 



provided r„ ^ 

Finally, in Fig. 4 the calculated temperature dependence 
of the electron-electron scattering time is depicted. Dotted 
curve represents a single- valley system, solid curve shows 
the calculation for the two-valley QW with zero valley- 
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orbit splitting and the same concentration of carriers. We 
remind that according to Eq. (15) the spin relaxation rate 
is obtained as a product of r* defined by Eq. (16) and the 
squared spin precession frequency taken at the Fermi level 
at zero temperature. The qualitative behavior of these two 
curves is similar: with the temperature increase the scat- 
tering time shortens as r* cx T"^ [see Eq. (21)] due to the 
weakening of Pauli blocking and reaches a minimum (seen 
in the figure only for the two-valley structure) caused by 
the transition to the non-degenerate case. This transition 
takes place at a smaller temperature for the two-valley 
system because the carrier concentration in each valley is 
twice smaller. For the accepted parameters the scatter- 
ing time in the two-valley system, in comparison with the 
single-valley system, is shorter at lower temperatures and 
longer at higher temperatures. 

One can sec from Fig. 4 that the scattering time r* 
has a picosecond scale in a wide range of temperatures. 
In the state-of-the-art Si/SiGe QWs where the spin re- 
laxation was studied the momentum scattering time Tp 
was about 10 ps for even smaller carrier concentrations 
than those taken in our calculation. Therefore, electron- 
electron collisions play a substantial role in controlling the 
spin relaxation in those Si/SiGe structures. 

Conclusions. — We have developed a theory of 
electron-electron scattering effect on the Dyakonov-Perel' 
spin relaxation in multi- valley semiconductor QWs. We 
have shown that, although the intervalley scattering of 
electrons is suppressed, the interaction of electrons occu- 
pying different valley-orbit split subbands influences the 
spin relaxation. The electron-electron scattering rates in 
single and multi-valley systems are different due to (i) re- 
distribution of electrons between the subbands, and (ii) 
an enhancement of screening in the two- valley systems. 

The values of electron-electron scattering times in high- 
mobility Si/SiGe QWs with two occupied valleys may be 
comparable and even shorter than the momentum scatter- 
ing time in a wide range of temperatures. Therefore, in 
these structures the electron spin relaxation can be con- 
trolled by electron-electron scattering. 
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